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present a Monte-Carlo AO simulation package that has been used to model 
the performace of EAGLE, and provide results, including comparisons with 



^ . an analytical code. These results include an investigation of the performance 
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O . cantly reduce the complexity of a real-time control system when implemented. 
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optics 

1. Introduction 

The next generation of optical ground-based Extremely Large Telescopes (ELTs) 
is currently in the design phase, with plans for primary mirror diameters of over 



30 m [1,2]. Once built, these facilities will allow astronomers to probe the universe 
with unprecedented sensitivity and very high resolution. A suite of instruments for 
these telescopes is planned, allowing many different observation goals to be met. 

The EAGLE instrument for the planned 42 m European ELT (E-ELT) is currently 
in the design phase [3]. It is a multi-object integral field unit spectrograph using 
adaptive optics (AO) with a multi-object AO (MOAO) system to correct incoming 
wavefronts in open-loop, using wavefront sensors which do not sense the corrections 
made to the science fields. Baseline designs for EAGLE include up to 11 wavefront 
sensors, using laser and natural guide stars. It is envisaged that there will be 20 science 
field pick-offs, allowing good AO correction in 20 separate fields, each 1.5 arcseconds 
diameter, simultaneously across a five arc-minute field. The use of a multi-object AO 
system allows good atmospheric correction to be achieved for selected objects across 
a wide field of view. 

Part of the design phase for EAGLE includes extensive simulation and modelling of 
the AO performance, since AO is an essential part of the instrument design. The sim- 
ulation and modelling are carried out in two phases. First, an analytical code is used 
to obtain an order-of-magnitude performance estimate, covering a large parameter 
space relatively quickly. However, many fine details which are essential to include for 
ELT scale designs are not included. A Monte-Carlo code is then used to fill in details 
giving a more reliable performance estimate, including non-linear effects, and noise 
sources. However Monte-Carlo simulation has far greater computational requirements, 



so a reduced parameter space is considered. 

The Durham AO simulation platform (DASP) is a Monte-Carlo code which can be 
used for the simulation of any common form of AO system (including classical AO, 
laser tomographic AO, multi-conjugate AO and MOAO [4]), and has been developed 
specifically with ELT simulation in mind [5]. It is an end-to-end time-domain code 
and is parallelised, allowing it to be used across a computing cluster using the Mes- 
sage Passing Interface (MPI) library to reduce computation time. It includes detailed 
models of telescope and AO systems, allowing high fidelity models to be produced. 

The development of a real-time control system for EAGLE is a challenge. There are 
expected to be of order 10^ wavefront slope measurements, and of order 10^ deformable 
mirror (DM) actuators to control per science path (of which there will be about 20). 
These wavefront slope measurements will be used to update the DM actuators at 
about 250 Hz. It is likely that EAGLE will use a conventional matrix-vector based 
wavefront reconstruction, though other techniques, such as iterative algorithms have 
not been ruled out. In this paper, we concentrate only on the matrix- vector based 
wavefront reconstruction and consider some details that may make this easier to 
implement in hardware. In order to access all the elements of the control matrix 
for each new set of slope measurements, a data rate of order 20 TBs~^ is therefore 
required (assuming four bytes per matrix element). This will require a very advanced 
control system, and so any achievable simplifications are desirable. 

Here, we present some recent results obtained from the simulation of EAGLE using 



the DASP. Some of these results are compared with those from an analytical code, 
produced independently by another member of the EAGLE consortium, where appro- 
priate. However, the analytical code is unable to include non-linear effects, and so can 
only be used for rough performance estimates. The simulations presented here include 
an investigation of compressed wavefront reconstructor algorithms which could sim- 
plify real-time control system design. The technical difficulty of the real-time control 
system design is such that it should be considered even at the early design phases. We 
also discuss the issue of Shack-Hartmann sensor non-linearities for open-loop systems. 
In §2, we describe the simulations that have been carried out, in §3 we give results, 
and conclusions are made in §4. 

2. Simulation description 

There are several possible designs for EAGLE, with different laser guide star (LGS) 
and natural guide star (NGS) requirements, based on trade-offs between cost, perfor- 
mance and sky coverage. Here, we concentrate on a design with nine LGSs equally 
spaced around a ring with a 7.3 arcmin diameter and a single NGS with 16 x 16 
sub-apertures used for low order mode correction (tip, tilt, focus and astigmatism) 
which LGS sensors are easily measure usefully. We assume that the LGSs are centre- 
launched and have an elongation of 5 arcsec at the edge of the telescope pupil (max- 
imum elongation). No measures to mitigate this elongation are made and wavefront 
slope computation uses a centre of gravity algorithm. The telescope diameter is as- 



sumed to be 42 m, and the wavefront sensors (WFSs) have 84 x 84 sub-apertures each 
with 20 x20 pixels, unless otherwise stated, requiring a WFS detector with 1680 x 1680 
pixels. Each science field pick-off has its own deforniable mirror (85 x 85 actuators) 
with wavefront control optimised along the line-of-sight for this field. Unless other- 
wise stated, the results presented here are for a target at the centre of the field, i.e. in 
the middle of the LGS ring. The science wavelength is H-band (1.65 /im). We use a 
virtual DM formulation for wavefront control, placing virtual DMs conjugate to the 
height of strong turbulence. This allows us to reconstruct the atmospheric turbulence 
at the positions of these virtual DMs, and use this knowledge to determine the shape 
that should be given to the physical DM for this science field. In the simulations here, 
the shape given to each MOAO physical DM is the sum of the virtual DMs projected 
along the line-of-sight for this science field. The virtual wavefront reconstruction is 
performed using a standard truncated least-squares matrix-vector algorithm with the 
vector containing the latest wavefront slope measurements, and the matrix being the 
pseudo-inverse of the system interaction matrix (the measured WFS response to per- 
turbations induced on the virtual DMs). The simulations presented here all assume 
two discrete layers of atmospheric turbulence, and two virtual DMs unless otherwise 
stated. It should be noted that the estimated performance reported here should be 
seen as optimistic due to the simple nature of a two layer profile. The atmospheric 
outer scale is taken as 50 m, and Fried's parameter is 10.6 cm (at 500 nm), corre- 
sponding to a seeing of 0.95 arcsec. These values are as used in simulations carried out 



by Fusco et. al. for the EAGLE consortium [6] which our Monte-Carlo simulations are 
used to verify. The update rate of the AO loop is 250 Hz, and unless otherwise stated, 
a delay (latency) of 4 ms between wavefront detection and correction is simulated. 

The deformable mirrors are operated in open-loop, i.e. the WFSs do not sense 
changes made to the DMs. 

In these simulations, atmospheric phase screens are translated across the telescope 
pupil assuming a frozen- flow turbulence model [7] . The sections of these screens rel- 
evant to a given line-of-sight at a given time are then selected (with sub-pixel in- 
terpolation) and summed (with interpolation for a source at finite distance, e.g. a 
LGS). These line-of-sight pupil phases are then used as input to Shack-Hartmann 
WFS models, which produce a simulated noisy Shack-Hartmann image, and to gen- 
erate science camera images, before and after correction of the phase using a DM. The 
DMs are controlled by a wavefront reconstructor, which uses the slope measurements 
taken by the WFSs to compute the correction to be applied. We use a centre of grav- 
ity algorithm for wavefront sensing. Sodium laser spots (as produced by a LGS) are 
assumed to form at 90 km with a depth of about 10 km and a Gaussian distribution. 
The performance of the AO system can be measured as a function of time, and the 
average long-exposure performance is also obtained. These simulations include many 
noise sources including detector noise, photon shot noise, laser guide star elongation, 
and WFS non-linearities. The simulation code is therefore suited to the high fidelity 
modelling of AO systems. 



2. A. Parameter space 

We have covered a large parameter space during these simulations, and with our 
available hardware we are able to cover about 7 parameter points per day including 
generation of interaction and control matrices. Parameters that have been explored 
include: 

1. Deformable mirror mis-conjugation 

2. Control matrix representation (investigating reductions in control matrix size 
to simplify real-time control system development) 

3. WFS linearisation 

4. Deformable mirror mis-alignment 

5. Wavefront sensor pixel scale 

6. Zenith angle 

7. LGS power 

8. WFS read-out noise 

9. Secondary mirror support obscuration 

10. Woofer- tweeter configuration 

11. Sodium layer profile 



Here, we consider further the first three parameters. DM mis-conjugation can occur 
when knowledge of the atmosphere is not perfect, and introduces an additional error 
into the correction applied to the DM. Control matrix representation is an important 
consideration for the design and development of real-time control systems, allowing 
designs to be simplified and costs reduced when the memory required to store a con- 
trol matrix is reduced. WFS linearisation is necessary for open-loop systems because 
Shack-Hartmann based WFSs have a slightly non-linear response to incident wave- 
front slope, and attempts to calibrate and correct this non-linearity can improve AO 
performance. 

3. Results 

3. A. Correction across the field of view 

Since MOAO systems operate in open-loop (the wavefront sensors do not sense the 
applied wavefront corrections), the corrections applied to the wavefront can be made 
along any line-of-sight. Unless stated otherwise, results presented here are for a line- 
of-sight at the centre of the field of view, i.e. the direction corresponding to the centre 
of the LGS ring. However, it is instructive to compare expected performance across 
the field of view, and Fig. 1 shows performance as a function of position across the 
field. In this figure, the corrected line-of-sight is moved from the on-axis location in 
a direction towards and past one of the LGSs (at 219 arcseconds). It can be seen 
that correction is uniform for most of the field of view within the LGS ring, and 



performance falls once the line-of-sight is close to the LGS ring, due to poor sampling 
of turbulence at these locations; turbulence here is only sampled by one WFS and 
so cannot be reconstructed well, while turbulence in the centre of the field of view is 
sampled by many WFSs. 

[Fig. 1 about here.] 

3.B. DM mis-conjugation 

A multi-conjugate AO system (including a MOAO system such as EAGLE, using 
virtual DMs) requires information about the strength and position of turbulent at- 
mospheric layers to operate most effectively. DMs are then conjugated at the locations 
of the most dominant layers. However, if mis-conjugation occurs, for example because 
layer positions are not well known, the AO system performance will be degraded. 

Fig. 2 shows how the performance of an AO system (correcting at H-band) is 
degraded by mis- conjugation. Here, dominant turbulent layers were placed at km 
and 10 km, and two virtual DMs placed at km, and at a varying height between 
8-12 km. This figure demonstrates that it is necessary to be able to conjugate DMs to 
within a few hundred meters of dominant turbulence. Analytical results provided by 
Fusco et al [6], which replace the LGSs with NGSs (the analytical code cannot model 
cone effect or spot elongation) for a system with 110 x 110 sub-apertures per wavefront 
sensor (long dashes) are shown to be slightly optimistic for perfect conjugation when 
compared with equivalent Monte-Carlo results (solid curve), with 110 x 110 sub- 



apertures (for comparison purposes) for each wavefront sensor. When a cone effect 
due to the laser spot being at a finite distance (meaning only a cone of turbulence 
is sampled by the WFS) and spot elongation caused by the three-dimensional nature 
of sodium emission are included in the simulations (dotted curve), performance is 
seen to fall by about ten percent. These effects are not modelled in the analytical 
results. A reduction of WFS order to 84 x 84 sub-apertures is shown to further reduce 
performance (the lower two curves). Performance is also shown to be dependent on 
the sodium profile. 

[Fig. 2 about here.] 

It should be noted that these simulations are a simphfication of the true situation 
where there will be many more turbulent layers, each with finite thickness. However, 
for all of these cases, the general trend with mis-conjugation is clear, implying that a 
DM should be conjugated to dominant turbulence with an accuracy of a few hundred 
meters. This places constraints on the design of turbulence profiling systems. 

Fig. 4 shows how the simulated AO system performance falls as a function of num- 
ber of atmospheric layers in these simulations. Here, we have not sought to optimise 
the wavefront reconstruction in any way, using a simple truncated least-squares wave- 
front reconstructor. A virtual DM has been placed conjugate to each layer, with an 
actuator spacing calculated to minimise fitting error, following [8]. The ideal number 
of virtual DMs, their conjugate heights and the actuator spacings to use to optimise 
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MOAO system performance is a subject of on-going research. Here, we do not con- 
sider the effect of DM mis-conjugation when there are more than two atmospheric 
layers. Table 3 shows the parameters used for these multiple layer simulations, as pro- 
vided by Fusco et. al. [6]. A global Fried parameter of 10.6 cm and an outer scale of 
50 m were used. It should be noted that wavefront reconstruction uses a least squares 
algorithm. The use of a minimum variance wavefront reconstruction may improve 
performance. However, this shows that the performance of EAGLE is likely to fall 
when the atmospheric turbulence is heavily layered. 

[Fig. 3 about here.] 

[Fig. 4 about here.] 

3.C. Reconstructor representation 

The control matrix for a single EAGLE science path is likely to contain of order 10^ 
elements, and must be accessed at a rate of 250 Hz, requiring a memory bandwidth 
of 1 TBs"^ assuming 32-bit floating point format storage. When considering that 
EAGLE is likely to have up to 20 science paths, the memory bandwidth requirement 
increases by a factor equal to the number of science paths, up to 20 TBs~^ for EAGLE. 
Reducing this memory bandwidth requirement is important to reduce the real- 
time control system complexity. Assuming a field programmable gate array (FPGA) 
based wavefront reconstruction unit, the memory bandwidth will be determined by 
the FPGA clock rate, the memory to FPGA bus width and the number of FPGAs 
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used for processing. By reducing the total size of the control matrix, the number of 
FPGAs can be reduced leading to a cheaper, simpler, more reliable design. We now 
consider several techniques that can be applied to reduce the control matrix size. 

3.C.I. Sparse representation 

Sparse matrix representation of AO system control matrices has been studied [9], 
and for multi-conjugate systems (or most systems without a specific WFS to DM 
alignment), sparse matrix techniques are known to perform poorly [10] due to poorly 
sensed modes and LGS tip-tilt uncertainty. Fig. 5 verifies this, showing that a highly 
non-sparse representation is required to maintain the AO system performance. The 
sparse matrices used here are created by removing the least infiuential parts of the 
control matrix, i.e. elements closest to zero. Typically, 70 % of the original matrix 
must be present, as demonstrated in Fig. 5. However, when stored in sparse format, 
each matrix element must be stored accompanied by its position, resulting in twice 
as much storage (assuming 32 bit fioating point for the matrix element, and a 32 bit 
integer for position), thus consuming more memory than the original control matrix. 
Therefore, sparse matrix representation is not a solution for EAGLE. 

[Fig. 5 about here.] 

3.C.2. Fixed point representation 

Fixed point representation is often used in hardware (for example FPGAs) as it is 
simpler to use than fioating point representations, and has a lower computational 
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complexity. Here, values are stored in twos-compliment integer format, with a known 
scaling factor. By using a fixed point control matrix representation, as shown in 
Fig. 6, it is possible to reduce the control matrix storage requirements by a factor 
of two, using 16 bit fixed point values rather than 32 bit floating point values, while 
still maintaining the AO system performance. To compute the fixed point control 
matrix, the minimum and maximum elements were found, and used to compute an 
offset (equal to the minimum value) and scaling factor (equal to the range), unique 
for a given control matrix. The fixed point control matrix elements are computed by 
subtracting the offset and dividing by the range before being scaled by 2^ where b is 
the number of bits used to store the fixed point representation. 

[Fig. 6 about here.] 

3.C.3. Compressed floating point format 

A control matrix is far from homogeneous, with a large range of values. This suggests 
that fixed point representation may not be ideal, since the relative resolution of small 
values is low and so will influence the wavefront error to a greater extent. We therefore 
consider a compressed floating point representation, which is able to cover the full 
range of 32 bit floating point, but with a reduced precision. Standard IEEE 32- 
bit floating point values have 8 bits dedicated to the exponent, 23 bits dedicated 
to the mantissa and a single sign bit. A compressed floating point format which 
reduces the precision of the mantissa can be investigated. In AO, wavefront slope 
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measurements are commonly computed using a centre of gravity measurement, which 
in good conditions (high hght level, low detector noise) is typically assumed accurate 
to at best a one hundredth of a pixel, and in practise, is far less accurate. With say 
20 X 20 pixels per sub-aperture, we can assume that there are 2000 measurable spot 
positions across the sub-aperture, which can be encoded in eleven bits. Therefore we 
can predict that a mantissa of a compressed floating point number need be no more 
than eleven bits wide. 

By running AO simulations with a range of bit-widths for the mantissa, we find 
(Fig. 7) that AO system performance is not degraded until fewer than 10-12 bits are 
used for the mantissa, which is represented by a compressed floating point number 
requiring between 19-21 bits in total. However, this is a greater storage requirement 
than we have shown to be required using a fixed point representation. 

[Fig. 7 about here.] 

3.C.4. Variable precision floating point format 

By compressing the exponent, as well as the mantissa of a floating point number, we 
can further reduce the storage requirement for the control matrix. We represent a 
floating point number in the form 

{-ly xbxa^ X (-+m) (1) 

where s is the sign (one bit), a is the base (2 for standard floating point representation) 
which is constant for a given control matrix, 6 is a scaling factor (constant for a given 
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control matrix), e is the exponent value, and m is the stored mantissa value. As with 
standard floating point representation, the mantissa is stored without an implicit 
integer part, which can be assumed (if it was not there, the exponent value can 
always be changed to shift the mantissa), and this is represented in the equation 
by the addition of the mantissa (flxed point with a value less than |) with |. The 
exponent, e is in standard twos-compliment integer format. 

To convert a standard control matrix into this format, the minimum and maximum 
values required for storage are flrst obtained. We then set requirements that the 
mantissa for the maximum value is all ones, and the mantissa for the minimum (non- 
zero) value is all zeros except for the flnal bit, which is set. The exponent for the 
maximum value has all bits set, and the exponent for the minimum value has all bits 
unset. A value of zero is represented by having all bits of the mantissa and exponent 
unset. These conditions allow us to flnd two unknown values, a and b which will allow 
us to store this control matrix with highest precision. We then proceed to convert the 
standard control matrix into the variable precision representation. 

As can be seen from Fig. 8, using a four bit exponent and mantissa is sufficient 
for good AO system performance, i.e. a total of nine bits per control matrix element 
(including a sign bit). Similarly, a five bit exponent and three bit mantissa, and a six 
bit exponent and two bit mantissa also provide similar performance (taking nine bits 
per control matrix element). The memory storage requirement has therefore been 
reduced by almost a factor of four. We have not investigated the effect of using a 
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greater number of atmospheric layers and virtual DMs, though at most, this will 
increase the number of bits required slightly. 

[Fig. 8 about here.] 

By using variable precision floating point for storage of the control matrix, the 
memory bandwidth requirement can be reduced by a factor of almost four, which 
will greatly simplify the design of a real-time control system for EAGLE. Only a 
quarter of the FPGAs used by an uncompressed system would be required, with 
simplifications also made by reducing the number of inter-FPGA connections, and an 
increased reliability due to a reduced number of components. 

3.C.5. Implementation in FPGA 

Implementation of variable precision floating point format in an FPGA is trivial: A 
24+4 _ 255 element look-up table can be used to translate the stored 9 bit control 
matrix values (using the mantissa and exponent for the index into the look-up table) 
into standard 32-bit floating point values to which the sign can then be inserted. A 
standard floating point multiplication routine can then be used during the matrix- 
vector multiplication. 

The Virtex-6 family of FPGAs is the latest offering from the company Xilinx, one 
of the major manufacturers of these devices. This range includes devices with up to 
1200 input/output pins, with 37 MB internal memory in the FPGA (similar to a 
central processing unit (CPU) cache), and a clock rate of up to 1.6 GHz. The design 
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of a real-time control system could be carried out using internal memory only. In this 
case, to store 20 control matrices (one for each line-of-sight) each 4 GB in size (1 
billion 32-bit floating point values), would require over 2000 FPGAs. Using variable 
precision floating point can reduce this requirement to just over 600 FPGAs, though 
this is still an undesirably large number. 

Alternatively, we can use external memory connected to the FPGA pins. With 
standard 32-bit floating point storage, and a 1024 bit wide memory bus (1024 FPGA 
pins connected to memory), we can access 32 values each memory read. The remaining 
pins are reserved for the address bus, and inter-FPGA communications. Assuming 
that memory can be accessed at the full FPGA clock rate (1.6 GHz) we will achieve a 
memory bandwidth of about 5 x 10^° values per second. The requirement for EAGLE 
is a minimum of 5 x 10^^ values per second, so 100 FPGAs would be required. If 
however, variable precision floating point storage is used, we could access 114 values 
each memory read (with a 1026 bit wide FPGA bus), equating to 1.8 x 10^^ values 
per second, requiring 28 FPGAs for EAGLE, a far more attractive proposition to 
develop. 

In practice, the memory bandwidth requirement may be increased to reduce the 
AO system latency. Here, we have assumed a latency of 4 ms, equal to the frame 
time (at 250 Hz). To achieve a latency of 1 ms, we would require an improvement in 
memory bandwidth by a factor of four, which in turn would require 112 FPGAs to 
meet this requirement. 
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3.D. Wavefront sensor calibration 

Shack-Hartmann based wavefront sensors are slightly non-linear due to the pixelated 
nature of the detector meaning that position information is lost: The measured slope 
is not proportional to the actual wavefront slope across the sub-apertures. For closed 
loop AO systems, this is not a problem since the degree of non-linearity is small and 
because the measured wavefront slopes are minimised by the DM, a linearity approx- 
imation works well. However, for typical open-loop systems, this is more problematic 
since large uncorrected wavefront slopes can be measured. Therefore, the corrected 
wavefront (unsensed) will have some additional error due to this non-linearity. This 
error is enough to lead to reduced performance of the AO system, and is present 
regardless of the slope measurement algorithm used if this algorithm is linear (e.g. 
centre of gravity, matched filter and correlation algorithms). However, a suitable cal- 
ibration of the WFS can be carried out, measuring the WFS estimated response to 
a set of known incident wavefront slopes (introduced by a flat mirror on a tip-tilt 
stage). During AO system operation, the uncalibrated measured wavefront slope can 
then be used to infer the true (calibrated) wavefront slope by interpolating from the 
calibration data. This calibrated measurement can then be used to perform a more 
accurate wavefront reconstruction. 

We have performed Monte-Carlo simulations using this technique for WFS cali- 
bration using a centre-of-gravity slope measurement algorithm, and have investigated 



the number of calibration steps required for good AO performance. These simula- 
tions are based around the aforementioned EAGLE simulations. We have used Shack- 
Hartmann sub-apertures with 20 x 20 pixels each, and a pixel scale of 0.8 arcsec per 
pixel at a wavelength of 589 nm and as before, Fried's parameter is 10.6 cm. The 
WFS calibration is performed over the entire sub-aperture field of view. This large 
field of view is due to the need to detect the elongated LGS spots, and due to the 
higher dynamic range of the open-loop WFS (spots are measured in open-loop, so are 
not necessarily close to the centre of sub-apertures, as is usually the case for a closed 
loop system). Fig. 9 demonstrates the degree of non-linearity in the simulated Shack- 
Hartmann wavefront sensor showing the deviation of measured wavefront slope from 
the true slope as the true spot position moves across the sub-aperture. A true (phys- 
ical) sensor would display even more non-linearity due to imperfect optics. Fig. 10 
shows the performance improvements achieved with increasing calibration accuracy 
(number of calibration steps), demonstrating that this linearisation calibration is an 
important part of open-loop AO system operation. We see that in this case, at least 50 
slope calibration measurements are required to achieve best performance, each step 
corresponding to a spot shift of less than half a pixel. By performing this calibration, 
the Strehl ratio (relative to uncalibrated performance) is increased by over 25%, and 
so the design of an open-loop real-time control system should therefore incorporate 
this calibration step. It should be noted that the optimal number of calibration steps 
is dependent on the WFS spot size on the detector so will vary with instrument and 
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atmospheric conditions. 

[Fig. 9 about here.] 
[Fig. 10 about here.] 

3.E. TMT comparisons 

The Thirty Meter Telescope (TMT) project also has plans for a multi-object spectro- 
graph with AO, IRMOS [11,12]. The results presented here show that the estimated 
performance of these systems (taking into account the many unknowns in the de- 
signs), both estimating 50-60% ensquared energy in 50 mas. It should be noted that 
the results presented in this paper have been for energy within 75 mas. When we use 
our simulation models to measure energy within 50 mas, this is typically about 1-2% 
lower than the energy within 75 mas. This serves to strengthen the assumption that 
modelling of AO systems can yield reliable performance estimates. 

4. Conclusion 

We have performed full end-to-end Monte-Carlo simulations of an AO system for EA- 
GLE. Investigations reported here show that the atmospheric turbulence profile must 
be well known, with the heights of turbulent layers known to within a few hundred 
metres. We have also reported on an investigation of compressed reconstructor repre- 
sentations and find that it is possible to reduce control matrix memory requirements 
by almost a factor of four in the cases investigated, significantly reducing the com- 
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plexity of an FPGA based real-time control system. An investigation into the effect 
that the non-linearity of Shack-Hartmann based wavefront sensors has on AO sys- 
tem performance has also been carried out, demonstrating that a linearity calibration 
should be included in an open-loop real-time control system to improve performance. 

Acknowledgements 

This work is funded by the STFC 

References 

1. J. Nelson and G. H. Sanders, "The status of the Thirty Meter Telescope project," 
in "Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Se- 
ries," , vol. 7012 of Society of Photo- Optical Instrumentation Engineers (SPIE) 
Conference Series (2008), vol. 7012 of Society of Photo- Optical Instrumentation 
Engineers (SPIE) Conference Series. 

2. J. Spyromilio, F. Comeron, S. D'Odorico, M. Kissler-Patig, and R. Gilmozzi, 
"Progress on the European Extremely Large Telescope," The Messenger 133, 
2-8 (2008). 

3. C. J. Evans, M. D. Lehnert, J. G. Cuby, S. L. Morris, A. M. Swinbank, W. D. 
Taylor, D. M. Alexander, N. P. F. Lorente, Y. Clenet, and T. Paumard, "Science 
Requirements for EAGLE for the E-ELT," in "Adaptive Optical Components II. 
Edited by Holly, Sandor ; James, Lawrence. Proceedings of SPIE, Volume 141, pp. 



21 



120-124," , vol. 7014 of Presented at the Society of Photo- Optical Instrumentation 
Engineers (SPIE) Conference (2008), vol. 7014 of Presented at the Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference, pp. 1-2. 

4. N. Devaney, "Review of astronomical adaptive optics systems and plans," in "So- 
ciety of Photo-Optical Instrumentation Engineers (SPIE) Conference Series," , 
vol. 6584 of Society of Photo-Optical Instrumentation Engineers (SPIE) Confer- 
ence Series (2007), vol. 6584 of Society of Photo- Optical Instrumentation Engi- 
neers (SPIE) Conference Series. 

5. A. G. Basden, T. Butterley, R. M. Myers, and R. W. Wilson, "Durham extremely 
large telescope adaptive optics simulation platform," Applied Optics 46, 1089- 
1098 (2007). 

6. T. Fusco, F. Assemat, E. Gendron, G. Rousset, N. A. Dipper, T. J. Morris, 
R. M. Myers, and B. Neichel, "EAGLE AO analysis and design report," EAGLE 
consortium document (2008). 

7. V. I. Tatarski, Wavefront Propagation in a Turbulent Medium (Dover, 1961). 

8. D. Gavel, R. Dekany, B. Bauman, and J. Nel- 
son, "Multi-conjugate deformable mirror fitting error," 
http://cfao.ucolick.org/research/aoforelt/MCAO_DM_FittingError.pdf (2003). 

9. W. J. Wild, E. J. Kibblewhite, and R. Vuilleumier, "Sparse matrix wave-front 
estimators for adaptive-optics systems for large ground-based telescopes," Optics 



22 



Letters 20, 955-957 (1995). 

10. B. L. Ellerbroek, "Efficient computation of minimum-variance wave-front recon- 
structors with sparse matrix techniques," J. Opt. Soc. Am. A pp. 1803-1816 
(2002). 

11. D. Gavel, B. Bauman, R. Dekany, M. Britton, and D. Andersen, "Adaptive op- 
tics designs for an infrared multi-object spectrograph on TMT," in "Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference Series," , vol. 6272 
of Society of Photo- Optical Instrumentation Engineers (SPIE) Conference Series 
(2006), vol. 6272 of Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series. 

12. D. R. Andersen, S. S. Eikenberry, M. Fletcher, W. Gardhouse, B. Leckie, J. Veran, 
D. Gavel, R. Clare, R. Guzman, L. Jolissaint, R. Julian, and W. Rambold, "The 
MOAO system of the IRMOS near-infrared multi-object spectrograph for TMT," 
in "Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Se- 
ries," , vol. 6269 of Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series (2006), vol. 6269 of Society of Photo- Optical Instrumentation 
Engineers (SPIE) Conference Series. 



23 



List of Figures 

1 A figure showing AO performance for different line-of-sigfit directions 
across tlie field of view. Tlie centre of tlie field of view is at arcsec, 

and an LGS is at 219 arcsec 25 

2 A figure showing the effect of DM mis-conjgation on AO system per- 
formance, which is represented by ensquared energy falling in a 75 mas 
diameter box with a science wavelength of 1.65 microns. The long- 
dashed curve shows analytical results for a 110 x 110 sub-aperture sys- 
tem, with the solid curve showing the Monte-Carlo simulation equiva- 
lent. The dotted curve just below this shows the performance reduction 
when a more realistic simulation including cone effect and LGS spot 
elongation is included. The two lowest curves show the performance 
reduction when the sub-aperture order is reduced to 84 x 84, with two 
different sodium layer profiles. misconjugationFl.eps 26 

3 A table showing the layer heights and relative strenghts used for mul- 
tiple layer simulations 27 

4 A figure showing AO system performance as a function of number of 
atmospheric layers and virtual DMs simulated 28 

5 A figure showing predicted EAGLE AO system performance as the 
sparsity of the control matrix is altered with the sparsity factor repre- 
senting the fraction of the original control matrix present. Uncertainties 

are about two percent in Strehl ratio. sparseF2.eps 29 

6 A figure showing predicted EAGLE AO system performance with a 
fixed point control matrix representation. Uncertainties are about two 
percent in Strehl ratio. fixedF3.eps 30 

7 A figure showing predicted EAGLE AO system performance with a 
compressed floating point control matrix representation. Uncertainties 

are about two percent in Strehl ratio. comfpF4.eps 31 

8 A figure showing predicted EAGLE AO system performance with a 
variable precision floating point control matrix representation. The key 
gives the number of bits used for the exponent for each curve. Uncer- 
tainties are about two percent in Strehl ratio. varfpF5.eps 32 

9 A figure showing the non-linearity of a Shack-Hartmann wavefront sen- 
sor, after subtraction of the linear response. calibrationF6.eps .... 33 

10 A figure showing predicted EAGLE AO system performance (Strehl ra- 
tio) as the number of calibration steps is increased. A 40% Strehl ratio 
is achieved with no calibration (0 steps, not shown on the logarithmic 
scale). linearstepsF7.eps 34 



24 



56 
55.5: 

55 
54.5 

54 
53.5 

53 
52.5 

52 
51.5 



-i h 



50 100 150 

Off-axis angle / arc-seconds 



Fig. 1. A figure showing AO performance for different line-of-siglit directions 
across tlie field of view. Tlie centre of tlie field of view is at arcsec, and an 
LGS is at 219 arcsec. 
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Fig. 2. A figure showing the effect of DM mis-conjgation on AO system perfor- 
mance, which is represented by ensquared energy falling in a 75 mas diameter 
box with a science wavelength of 1.65 microns. The long-dashed curve shows 
analytical results for a 110 x 110 sub-aperture system, with the solid curve 
showing the Monte-Carlo simulation equivalent. The dotted curve just below 
this shows the performance reduction when a more realistic simulation includ- 
ing cone effect and LGS spot elongation is included. The two lowest curves 
show the performance reduction when the sub-aperture order is reduced to 
84 X 84, with two different sodium layer profiles. misconjugationFl.eps 
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Number of layer Layer heights Layer strengths 

2 0, 12800 0.92, 0.08 

3 0, 1800, 12800 0.77, 0.17, 0.06 

4 0, 1800, 4500, 12800 0.67, 0.15, 0.13, 0.05 

5 0, 300, 1800, 4500, 12800 0.53, 0.21, 0.12, 0.10, 0.045 

6 0, 300, 900, 1800, 4500, 12800 0.47, 0.18, 0.11, 0.1, 0.09, 0.04 



Fig. 3. A table showing the layer heights and relative strenghts used for mul- 
tiple layer simulations 
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Fig. 4. A figure showing AO system performance as a function of number of 
atmospheric layers and virtual DMs simulated. 
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Fig. 5. A figure showing predicted EAGLE AO system performance as tlie 
sparsity of the control matrix is altered with the sparsity factor representing 
the fraction of the original control matrix present. Uncertainties are about two 
percent in Strehl ratio. sparseF2.eps 
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Fig. 6. A figure showing predicted EAGLE AO system performance witli a 
fixed point control matrix representation. Uncertainties are about two percent 
in Strehl ratio. fixedF3.eps 
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Fig. 7. A figure showing predicted EAGLE AO system performance witli a com- 
pressed floating point control matrix representation. Uncertainties are about 
two percent in Strehl ratio. comfpF4.eps 
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Fig. 8. A figure showing predicted EAGLE AO system performance witli a 
variable precision floating point control matrix representation. The key gives 
the number of bits used for the exponent for each curve. Uncertainties are 
about two percent in Strehl ratio. varfpF5.eps 
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Fig. 9. A figure showing the non-hnearity of a Shack-Hartmann wavefront 
sensor, after subtraction of the linear response. calibrationF6.eps 
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Fig. 10. A figure showing predicted EAGLE AO system performance (Strelil 
ratio) as tlie number of calibration steps is increased. A 40% Strehl ratio is 
achieved with no cahbration (0 steps, not shown on the logarithmic scale). 
linearstepsF7.eps 
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